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ABSTRACT 

Far-infrared [C ii] 158-/ini, [O i] 145-/im and [O i] 63-/ini fine structure emission 
fine fluxes were measured from archival Infrared Space Observatory Long Wavelength 
Spectrometer spectra of 46 extragalactic sources, with 28 sources providing detections 
in all three lines. For 12 of the sources, the contribution to the [C ll] 158-/J,m line flux 
from H II regions could be estimated from their detected [N ii] 122-fim line fluxes. The 
measured [C ii]/[0 i] and [O i] 63/145-/im line flux ratios were compared with those 
from a grid of PDR models previously computed using the UCL_PDR code. Persistent 
offsets between the observed and modelled line ratios could be partly attributed to the 
effects of [O i] 63-yLtm self-absorption. Using the SMMOL code, we calculated model 
[O i] line profiles and found that the strength of the [O i] Q3-fim line was reduced 
by 20-80%, depending on the PDR parameters. We conclude that high PDR densities 
and radiation field strengths, coupled with the effects of [O i] 63-fim self-absorption, 
are likely to provide the best match to the observed line flux ratios. 

Key words: FIR fine structure emission lines, PDRs, Extragalactic Sources, ISO 
LWS, Oxygen self- absorption. 
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1 INTRODUCTION 

Studies of our own Galaxy and observations of external 
galaxies have suggested that stellar ultraviolet radiation can 
ionize vast volumes of a galaxy and that far-ultraviolet radi- 
ation impinging on neutral cloud surfaces is responsible for 
a large fraction of the observed far-infrared (FIR) sp ectral 
line emission that cools the gas (jCrawford et al. 

mm- The 

relative contributions from different gas phases can be esti- 
mated from observations of several FIR forbidden lines. Fine 
structure (FS) emission lines can be used as tracers of nebu- 
lar conditions such as density, excitation and ionization. By 
virtue of their different excitation potentials and critical den- 
sities, FS emission lines provide an insight into the energet- 
ics and chemical composition of the regions from which they 
originate. The FIR fine structure lines of abundant species, 
such as [C ii] and [O i], have long been recognised as one 
of the most important coolants in the interstellar medium 
(ISM). In particular, the emission from singly ionized car- 
bon [C ii], at A=158 /xm, is known to trace warm and dense 
photodissociation regions (PDRs) and is used to trace the 
gas surrounding active star-forming regions. Carbon is the 
fourth most abundant element and has a lower ionization 
potential (11.26 eV) than hydrogen. For this reason it is 
predominantly in the form of C"'' in the surface layers of 



PDRs, where photoionization of neutral carbon is effective, 
but can potentially also be excited in ionized regions. [C ll] 
158/im emission is reasonably easy to excite, with a tran- 
sition energy of 92 K. The line is widely observed and is 
usually optically thin. The depth of the zones in neutral 
clouds is generally determined by dust extinction of carbon 
ionizing photons and often extend s to At, ^4, alth ough C^ 
decreases in abundance for At,^l (|Bell et al.ll2006l ). 

Another important species is oxygen, which has an ion- 
ization potential of 13.62 eV, slightly above that of hydro- 
gen. This means that the ionization structure of 0+\0'' 
closely follows the ionization structure of hydrogen. Also, 
the charge exchange of neutral oxygen in two-body recom- 
bination with hydrogen cannot be neglected as it affects 
the ionization balance. The incident FUV radiation main- 
tains a significant abundance of atomic oxygen deep into 
the cloud through photodestruction of CO and O2. All oxy- 
gen not incorporated into CO can remain in atomic form 
to dept hs as large as A„ = 10 mag un der strong FUV irra- 
diation (ISternberg fc Dalgarnol ll995^ .The fraction of oxy- 
gen at magnitudes greater than 10 depends on the physi- 
cal conditions in the cloud (i.e. its density, metallicities and 
radiation field strength), but essentially if the gas density 
is >10^ cm~^ at these depths, the oxygen not locked in 
CO is depleted on the grains. However for the purpose of 
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these calculations, the code that we use, UCL_PDR, does 
not take into account freeze out reactions. Atomic oxygen 
has two fine structure transitions, at 63 and 145 /im. The 
[O l] line at 63/im has an excitation energy corresponding 
to 228 K, is emitted predominantly by warm and dense 
gas, although it can also h ave an ionized gas component 
l|Aannestad fc Emervl |2003| ). and can act as a coolant in 
neutral hydrogen regions. It may also be a heating agent 
as O I can heat the gas through the coUisional de-excitation 
of th e dust continuum radiation- excited ^Pi fine structure 
level iTielens fc Hollenbach|[l985l ) . However this mechanism 
is not well studied . This line is optically thick in some cases 
(|Abel et all |2007| ) and its optical depth can be estimated 
by comparing the [Q l] line flux at 63/im to that of the [O 
l] line at 145^im (jAbel et al.ll2007l l. The [O l] 145/im line 
can be harder to detect: in fact the relative faintness of this 
line diminishes its capabilities as a robust PDR diagnostic. 
The lower level of the 145/xm line is not the ground state 
of 0°, me aning that [0 l] 145Atm emission is usually opti- 
cally thin (|Abel et al.ll2007l ). The 145/im line has an excita- 
tion energy of 326 K, almost 100 K higher than that of the 
63pim line, so the ratio of the two line intensities is sensi- 
tive to the gas temperature for T^;300 K. Since the strength 
of the FUV radiation field governs the heating of the gas, 
this line ratio is also a diagn ostic of the incident FUV flux 
l|Tielens fc Hollenbach|[l985l ). 

In the presence of foreground cold and tenuous gas 
the [O l] 63/^m li ne can show both emission and absorp- 
tion c omponents (iTielens fc HollenbaciilllQSSl : |l iseau et al.l 
I2OO6I : iGonzalez- Alfonso et al.ll2004h . Nevertheless observa- 
tions of the FIR fine structure emission lines of [O l] and 
[C 11] can be used as diagnostics to infer the physical con- 
ditions in the emitting gas, such as temperatures, densities 
and radiation field strengt hs, by comparing with mod els of 
photodissoci ation regions (iTielens fc Hollenbachlll985^ and 
Hll regions (|Rubin et al.lll991 ). In the past, several models 
have been used to interpret infrared, submillimeter and mil- 
limeter line observations of neutral gas in our own Galaxy 
as w ell as in the centr a l regions of nearby starb u rst galaxies 
(e.g. iMao et all I2OO0I : ISeaauist fc Fraveil I2OO0I : IWUd et al.l 
I1992I . and references therein). 

In this paper we will use FIR line fluxes measured from 
the ISO archive for a sample of extragalactic sources. We 
present the fluxes of fine structure emission lines observed 
by ISO from 46 sources and we estimate the bulk properties 
of the ga s in these e xtern al galaxies using the UCL_PDR 
model of iBell et all (|2005|) and the radiative transfer code 
SMMOL (|Rawlings fc Yatesll200lh . 



The paper is structured as follows: in Section 2 we give 
details of the sources selected and the fine structure emis- 
sion line fluxes we measured. In Section 3 we describe the 
UCL_PDR models. In Section 4 we evaluate the sensitiv- 
ity of the UCL_PDR models to variations in the physical 
parameters. In Section 5 we will consider the contribution 
of H II regions to fine structure emission lines. The results 
are given in Section 6, and there we discuss their interpreta- 
tion in the context of our SMMOL line profile calculations. 
Concluding remarks are given in Section 7. 



Source Index Number 

Figure 1. [C n]i58 /[O l]63 and [O l]63/[0 l]i45 ratios. Filled tri- 
angles represent the ratio of the literature value of [O l]63/[0 l]i45 
over our own measured [O l]63/[0 l]i45 ratio, with respective er- 
ror bars. Open triangles represent the ratio of the literature value 
of [C n]i58 /[O l]63 over our own measured [C Il]i58 /[O l]63 ra- 
tio with respective error bars. The horizontal line represents the 
case in which literature ratios and our own measured ratios are 
perfectly matched. The source index numbers and galaxy identi- 
fications are listed in Table [T] 



2 THE ISO DATA 

All the FIR data included in this paper were ac quired using 
the I SO Long Wavelength Spectrometer (LWS; IClegg et al.l 
Il996h . which had an estimated FWHM beam size of ~ 80 
arcsec. 

We collected [C ll]158pim, [O l]145/im and 63/xm emis- 
sion line fluxes from ISO LWS archival observations for 46 
external galaxies. The [C ll]158/xm, [O l]145/im, 63/im emis- 
sion line fluxes were measured from the spectra using the 
interactive package SM ART (Spectroscopy M odeling Anal- 
ysis and reduction Tool; iHigdon et al.l (|2004f ).'l 

Table [1] lists the emission line fluxes measured by us 
from the ISO spectra. For 33 of the galaxies, fluxes pre- 
viously published in the literature were available and these 
are also listed in Table[T] When uncertainties were not listed 
for these literature fluxes, we have estimated them as 20% 
for fluxes 10"^"* W m"^ and 10% for fluxes ^ 10"" W 
vcC^ . Although the choice of how to estimate the associ- 
ated errors, when it is not provided, is arbitrary, we used 
typical uncertainties resulting from calibration and pointing 
error s in the submillimeter and IR dom ain (|Martfn et al.l 
I2OO6I : iBavet et al]|2004 llsrael et al.lll995l ). The percentage 
uncertainty estimated in this way appears, on average, to be 
similar to the percentage line flux uncertainties measured by 
us using the interactive package SMART. 

In order to show that our measurements are in good 
agreement with the fluxes collected from the literature we 
have plotted our measurements versus literature values in 
Fig [1] We have plotted the ratio of the literature value of 
[C Il]i5g /[O l]63 over our own measured [C Il]i58 /[O l]63 
value, with respective error bars (open triangles), as well 
as the ratio of the literature value of [O l]63/[0 l]i45 over 
our own value of [O l]63/[0 l]i45, with respective error bars 
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Table 1. Line flux measurements, in units of 10 1* W m ^, for extragalactic sources with ISO-LWS detections of all 3 FIR fine structure 
emission lines of [C ll] and [O l] . For each source, our own line flux measurements are listed in the first row, while literature values, if 
any, are listed in the second row llContursi et al.ll20 02': 'Bcr gvall et al.ll20qol: IColbert et ailllQQgl : iFischer et aLlll996l : iLuhman et~al]|2003l ; 
iMalhotra et al.ll200ll : iNegishi et al.ll200ll : lUnger et al.„2000i : iBrauher et al.ll2008l'l . (n.a: not available) 
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<0.016 


1.41ib0.26 


n.a 










0, 


.40±0.08 


0.18±0.04 


n.a 


2.22ib0.93 


n.a 


34 


NGC 6946 


48 


45700139 


0, 


.876±0.028 


0.611±0.024 


<0.061 


1.43ib0.10 


n.a 










1 


.03±0.21 


0.59±0.12 


O.OSiO.Ol 


1.74ib0.71 


11.8±4.8 
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Index 


Source 


Radial Vel 

(km s^l) 


TDT No. of 
Observation 


[CII] 158fim 


[OI] 63(jni 


[OI] 145fim 


[CII] 158 
[OI]63 


[OI]63 
[OI]i45 


35 


NGC 


7673 


3408 


76601364 


0.1061±0.0041 


0.172±0.041 


<0.0075 


0.62±0.17 


n.a 












n.a 


n.a 


n.a 


n.a 


n.a 


36 


Mrk 331 


5541 


56500637 


0.182±0.006 


0.096±0.008 


0.0076±0.0018 


1.89±0.22 


12.63±4.03 












0.148±0.008 


0.11±0.006 


< 0.13 


1.34±0.14 


n.a 


37 


NGC 


4151 


995 


35300163 


0.822±0.006 


0.445±0.024 


<0.042 


1.85±0.11 


n.a 












0.074±0.006 


0.37±0.029 


n.a 


0.20±0.03 


n.a 


38 


NGC 


6286 


5501 


20700509 


0.187±0.008 


0.087±0.004 


0.0094±0.0013 


2.15±0.19 


9.25±1.70 












0.168±0.005 


0.073±0.009 


0.005±0.001 


2.30±0.35 


14.6±4.7 


39 


NGC 


6574 


2282 


70500604 


0.436±0.021 


0.223±0.034 


0.078±0.017 


1.95±0.39 


2.86±1.05 












0.44±0.020 


0.23±0.041 


n.a 


1.91±0.42 


n.a 


40 


NGC 


6822 


-57 


34300915 


0.176±0.011 


0.128±0.005 


0.008±0.002 


1.37±0.15 


16.5±3.9 












0.187±0.011 


0.143±0.006 


<0.05 


1.30±0.13 


n.a 


41 


NGC 


7552 


1608 


36903087 


0.583±0.028 


0.20±0.04 


0.037±0.007 


2.92±0.72 


5.41±2.10 












0.64±0.015 


0.63±0.023 


n.a 


1.01±0.06 


n.a 


42 


NGC 


7771 


4277 


56500772 


0.29±0.01 


0.141±0.009 


0.0187±0.0026 


2.06±0.20 


7.54±1.52 












0.298±0.009 


0.115±0.013 


n.a 


2.59±0.37 


n.a 


43 


NGC 


4041 


1234 


22202506 


0.32±0.01 


0.38±0.05 


<0.14 


0.84±0.13 


n.a 












0.348±0.005 


0.197±0.016 


n.a 


1.76±0.16 


n.a 


44 


NGC 


0278 


627 


59702260 


0.72±0.002 


<0.03 


0.020±0.002 


n.a 


n.a 












0.728±0.012 


0.347±0.016 


<0.021 


2.09±0.13 


n.a 


45 


NGC 


0695 


9735 


63300744 


0.137±0.007 


0.19±0.01 


0.010±0.002 


0.72±0.07 


19±4.8 












0.20±0.008 


0.118±0.009 


n.a 


1.69±0.19 


n.a 


46 


NGC 


0986 


1974 


74300187 


0.278±0.009 


0.19±0.04 


0.013±0.001 


1.46±0.35 


14.62±4.20 












0.304±0.011 


0.127±0.013 


<0.025 


2.39±0.33 


n.a 





1 1 1 1 1 1 1 


1 1 1 1 1 1 





















-1.0 -0.5 0.0 0.5 1.0 

Log(CII/OI63) 

Figure 2. The observed [C Il]i58 /[O l]63 ratios versus the 
[O l]63/[0 l]i45 ratios for detections ^ 4(t, with error bars. 



(filled triangles). The horizontal line in Fig [T] represents the 
case in which the literature ratios and our own measured 
ratios are in perfect agreement. 

Flux ratios of two lines of the same species, such as 
[O l]i45/[0 l]63 can provide information about nebular con- 
ditions such as temperature (or radiation field strength). In 
addition, a ratio such as [O in]gg/[0 ill] 52 can provide infor- 
mation on the electron density. In addition by using ratios, 
rather than fitting the line intensities directly, the beam fill- 
ing factors of the two emission lines cancel out, assuming 
that they come from the same regions and are the same 
size. 

Fig [2] shows the [C 11] 158 /[O l]63 ratio versus the 
[O l] 63/(0 l]i45 ratio for the 28 sources, out of the 46 listed 
in Table [1] that showed 4a" detections in all three lines. 



3 THE UCL_PDR MODELS 



We us ed the UCL_PD R time and depth-dependent PDR 
code (|Bell et al.l I2OO6I ) which solves the chemistry, ther- 
mal balance and radiative transfer self-consistently within 
a cloud. The range of parameter values input to the code 
are listed in Table [S] including the mean grain radius and 
grain number density adopted for grain surface reactions. 
The FUV g rain opacity prop erties adopted by UCL_PDR 
are listed bv lBell et af] (|2006l l. The outputs of the code are 
the abundances of each species included in the chemical net- 
work, their column densities, the emissivities and integrated 
intensities of the emission lines involved in the cooling, and 
the gas and dust temperatures. All of these quantities are 
functions of depth and time. At each depth-step, the model 
calculates the attenuation of the FUV field before begin- 
ning an iterative cycle to determine the gas temperature at 
which the total heating and cooling rates are equal, that is, 
when the condition of thermal balance is satisfied. For each 
iteration, the chemistry is first calculated, based on the gas 
temperature and attenuated FUV flux, after which the heat- 
ing and cooling rates are computed, together with the ra- 
diative transfer in the emission line, using the revised chem- 
ical abundances. The difference between the total heating 
and cooling rates is then used to derive a new gas tempera- 
ture. We use a grid of 1702 PDR models, already computed 
and partially used in previous works (Bell at al. 2006, Bell 
& Hartquist et al. 2006, Bell et al. 2007) spanning a large 
range of densities, metallicities and cosmic ray ionization 
rates. Amongst the range of parameter values covered by 
the grid of models (see Table [2]) , we do not consider vari- 
ations of the turbulent velocity parameter in our attempts 
to fit the observed emission from each galaxy, since the line 
ratios are believed to be fairly insensitive to small changes 
in the turbulent velocity jWolfire et al.lll989l l. 
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Table 2. Physical and chemical parameters for the UCL_PDR 
grid. In the grid of models the density and the radiation field 
parameters are incremented by 2 dex and 1 dex, respectively. 



Parameter Range of values Reference 

parameter 
values 

Cloud density (cm 10^^ nn $10^ nH= 10^ 

Incident FUV flux (Habing) 10^ Gq^IO^ Go = 10 

Age of the cloud (yr) lO^^tsJlO^ t = 10'' 

Cloud size (mag) 0^ A^,^10 A^ = 10 

Mctallicity Z/Zq =5,4,3,2, 1,0.5 Z=Zq 

0.1,0.25,0.01 

C.R. ionization rate (s~^) C = 5x,50x,500xlO^^'' C=5xia^^'' 

5x,50x,500xl0^^'= 
5x,50x,500xl0~^'' 

Turbulent velocity (km s ^) — — 

Carbon elemental abundance 1.4x10 ^n/f Sofia Sz Myer 

(2001) 

Oxygen elemental abundance 3.2xl0~^n/f Sofia Sz Myer 

(2001) 

Grain radius 0.1 0.1 ^m 

Grain number density Ug 2xlO~^^Z n cm^'^ 2x10"^^ cm"""^ 

4 SENSITIVITY OF MODEL RESULTS TO 
PARAMETER VARIATIONS 

We attempted to find the range of physical parameters wliicli 
best reproduced tiie observations, such as visual extinction, 
metallicity, cosmic ray ionization rate and UV radiation field 
strength. The influence of changing the number density of 
hydrogen nuclei (n_r/), the visual extinction (A„), the inci- 
dent radiation field strength (Go), the cosmic ray ionization 
rate ( and the metallicity (Z) is considered individually by 
varying only one of these parameters at a time. The age of 
the modeled regions was set to a value of 10^ yrs, owing 
to the fact that clouds with ages ^10^ yrs do not undergo 
significant changes in their predicted chemical profiles, even 
though, in some cases, the c hemistry only re aches its final 
steady state somewhat later (|Bell et al.ll2006l ). 

When one parameter was varied the remainder were 
held constant at the reference parameter values listed in the 
final column of Table [21 

In Figs 3 through 8 we overplot our model predictions 
against flux ratios from the [C ll] 158/im and [O l] 63/im 
and 145/im ISO LWS fluxes. Table [3] lists these 3 important 
PDR cooling transitions, their wavelengths, upper energy 
levels Eupper, and critical densities ricr- In this section we 
will only discuss the sensitivity of the ratios to changes in the 
input parameters. The comparison with observations will be 
made in Section 6. 



4-0.1 Visual Extinction 

PDRs derive their properties primarily from the penetration 
of FUV photons into their interiors. Dust grains provide the 
major source of continuum opacity and attenuate the inci- 
dent radiation field by selectively absorbing and scattering 
light at visible and ultraviolet wavelengths. This means that 
there is a strong correlation between visual extinction and 



Table 3. PDR diagnostic transitions. Th e critical d ensities for 
[C ll] and [O l] are for collisions with H l|TieIens fc HollcnbachI 
Il985l) . 



Species 


Transitions 


Wavelength 
(Mm) 


(K) 


(cm-3) 


[C II] 


^P3/2 Pl/2 


157.74 


92 


3x10^ 


[O I] 


3p 3 p 


63.18 


228 


4.7x10^ 


[O I] 


3P0 Pi 


145.53 


326 


1x10^ 



the chemistry of PDRs. We examined models for different 
visual extinction values in the range 1^A„^10. From Fig[3] 
one can see that there is no significant change in the emission 
line ratios from Ai, = l (black circle) to Ai, = 10 (green circle). 
This trend implies that the majority of the line fiuxes arise 
from 1^A„^2 and is consistent with the C"*" abundance be- 
cause, for values of A^ larger than 2, this ion recombines to 
form C" and subsequently forms CO. Indeed is a tracer 
of the edge of PDRs. This implies that the C"^ fine struc- 
ture emission lines that are observed are coming predom- 
inantly from regions with low visual extinction. Although 
neutral oxygen is still present at higher visual extinctions, 
it is somewhat reduced beyond A„ = l b ecause a fraction o f 
oxygen combines with C to form CO (see iRolhg et al.ll2007l ). 
We fixed Ai, = 10 as the most appropriate depth range to fit 
our sample of sources, since Fig [3] shows that we cannot dis- 
tinguish between models with optical depths in the range of 
IsjA^^lO. 

4.0.2 Radiation field 

Radiation field strengths incident upon a gas cloud can 
vary from the stan dard interste llar field ( x~ Draine) up 
to x~10^ Drai nes (lDrainelll978l) in region s of intense star 
formation (see iHoUenbach fc Tielend [l997l ) . The FUV field 
may be also expressed in terms of the Habing parameter 
(|Habinelll968h . The FUV fiux expressed in this way is de- 
termined by the parameter Go^x/1-'''- We adopt a range of 
lO^Go ^500000 Habings for consideration in this study, but 
we only plot the significant results. From Fig|4l [C Il]i58 /[O 
i]63 versus [O l]63/[0 l]i45 appears to be sensitive to varia- 
tions in the radiation field strength. As the radiation field 
strength increases so does the [O l]63/[0 l]i45 ratio, while 
the [C 11] 158 /[O l]63 ratio decreases. The energy required to 
excite the [O l]63 transition is somewhat higher than that 
for [C 11] 158, hence the line intensity ratio [O l]63/[C 11] 153 is 
expected to increase with gas temperature and hence higher 
FUV flux. PDR models predict that for nH>W^ cm'^ the 
line ratio will increase with both Go and uh, due to the dif- 
feren t critical densities of the two transitions (jWolfire et al.l 
|1990D . 

4.0.3 Metallicity 

Metallicities significantly below solar are observed in Lo- 
cal Group galaxies, including the Small Magellanic Cloud, 
as well as in more distant dwarf galaxies, with I Zw 
18's metallicity of 1/40 being amongst the lowest known 
(|lzotov fc Thuar] Il999l ). The metallicity-dependence ap- 
pears in several key processes in the UCL_PDR code. El- 
emental abundances of all metals are assumed to scale lin- 
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5.6074x10"' 5,6076x10 



0.0 0.5 
Log(CII/OI63) 




Log(CII/0n45) 



Figure 3. [C lljisg /[O l]63 versus [O l]63/[0 l]i45 ratios (left) and [C n]i58 /[O l]i45 versus [O l]i45/[0 ijea ratios (right), with error 
bars, compared with models of varying visual extinctions. The insert, in the left panel, shows a zoomed-in region of the overlapping 
coloured circles corresponding to differing visual extinction values. 



early with metallicity (Z/Z©); the dust-to-gas mass ratio is 
also assumed to scale linearly with metallicity and takes a 
standard a value of 10"'^ at Solar metallicity. The forma- 
tion rate of H2 on grain surfaces and the grain photoelectric 
heating rate are assumed to scale linearly with metallicity. 

We investigated a metallicity range of 0.01^ Z/Zq^S 
for this study, with the adopted solar neighbourhood car- 
bon and oxygen elemental abundances listed in Table [S] 
Metallicity affects the total abundances of carbon and oxy- 
gen bearing species, and hence can influence the chemical 
and thermal structure of PDRs in galaxies. 

In Fig [5] we show our model results for metallicities 
of 0.01 (hght green circle), 0.25 (claret circle) and 0.5 
(turquoise circle) , with remaining symbols representing val- 
ues of 1, 2, 3, 4 and 5 times solar metallicity respectively. 

The C/O ratios for the Sun, the Large Magellanic Cloud 
(LMC) and the SmaU MageUanic Cloud (SMC) are all --0.5 
within the uncertainties, but for lower metallicities than the 
SMC C/H decreases faster than O/H. We tried to reproduce 
this trend by changing carbon and oxygen elemental abun- 
dances with metallicity. In Fig|6]we show ou r model results 
for a metallicity of Z/Zq=1, with C/O=0.53 (|Asplund et al.l 
l2005h . then we scaled all the elemental abundances with 
Z/Z0 down to Z/Z0=O.25. Below Z/Zq=0.25 the C/O ra- 
tio was scaled from 0.53 to 0.16 as Z /Zp) was reduced from 
0.25 to 0.03 (jlzotov fc ThuanI Il999^ . As expected, the [O 
i]63/[0 l]i45 ratio is not sensitive to variations of metallic- 
ity whereas the [C 11] 158 /[O l]63 ratio would appear at first 
sight to be a good tracer of metallicity. 



4.0.4 Density 

We have considered the sensitivity of the fine structure emis- 
sion line intensities to H-nuclei densities in the range 10^< 
uh^JIO^ cm^'^. In Fig[7]we show our model results for dif- 
ferent values of density. The [C Il]i58 /[O l]m and [O l]63/[0 
l]i45 ratios are both sensitive to density, however there is a 
higher sensitivity for the [C Il]i58 /[O i]63 ratio. As expected. 



there is a substantial decline in the [O l]63/[0 i]i45 ratio 
as the density decreases from nH-^lO^ cm""^ to nH~9xl0^ 
cm^"^, and a steady decline in the [O l]63/[0 l]i45 ratio as 
the density decreases from n_r/~9xl0^ while the model [C 
Ii]i5s /[O i]63 ratios increase initially with decreasing den- 
sity until the density reaches a value of uh^IO^ cm""^ and 
then decrease as the density decreases further. Clearly, the 
changes in the [O l]63/[0 l]i45 ratio at 10^ cm~^ are partly 
due to the critical density of the two lines (see Table 
Similarly, the decrease in the [C 11] 158 /[O i]63 ratio as the 
density decreases is partially due to the density falling below 
the critical density of the [C 11] 158 line. 



4.0.5 Cosmic ray ionization rate 

The cosmic ray flux is known t o vary by over an o rder of 
magnitude in the Milky Way (|Schiike et all 1 19931 ). Mag- 
netic field lines can channel cosmic rays away from dense 
molecular cores; alternatively, the fiux of particles in star 
forming re gions can be many times higher than the canon- 
ical rate (jschilke et all 1 19931 ). The UCL_PDR code does 
not have a separate treatment for X-ray Dominated Region 
(XDR) effects that m ay be important in some galaxies (see 
iMeiierink et al.ll2006l ). However to a first approximation one 
may use an enhanced cosmic ray ionization rate to mimic the 
effects of XDRs (BeU & Hartquist et al. 2006). In Fig E] we 
show our model results for cosmic ray ionization rates in the 
range 5xlO~^^^(^^55xlO~^^s~^. The higher en d of our range 
can r epresent a higher ionization rate (see e.g. iMcCall et al.l 
l2003h or the effect of additional ionization due to X-rays. 
There is a substantial decline in the [O l]63/[0 l]i45 ratio as 
the cosmic ray ionization rate decreases from (^=5xl0~^^s~^ 
(green circle) to ("=5xl0~^^s~^ (claret circle) and it then de- 
creases slightly for <^=5xl0^^^s^^ (blue circle); the [C Il]i58 
/[O i]63 ratio decreases significantly from (^=5xl0~^^s~^ 
(green circle) to (^=5xl0~^''s^^ (claret circle) and then re- 
mains constant. In fact both the [C 11] and [O l]63 fluxes 
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Figure 4. The four panels show observed [C Il]i58 /[O l]63 versus [O l]63/[0 l]i45 ratios, with uncertainties. In the top left panel the 
coloured circles represent models with radiation field strengths in the range 10'':CGo^^2xlO'*. In the top right panel the coloured circles 
represent models with radiation field strengths in the range lO^SjGo^SxlO"^. In the bottom left panel the coloured circles represent 
models with radiation field strengths in the range lO^^Go^SxlO^. In the bottom right panel the coloured circles represent models with 
radiation field strengths in the range lOSJCo^JSO. 



decrease with a decrease in ionization rate, with [O ijea de- 
creasing at a faster rate. 



5 USING [N II] FLUXES TO ESTIMATE THE 
IONIZED GAS CONTRIBUTION TO [C II] 
FLUXES 

Several studies have been made that show that significant 
[C ll] emission can arise from H ll regions (jHcilcs 1993 : 
lAbel et al.|[2005l : iKaufman et"aLll2006l : lAbellbOod ). Since the 
ionization potential from to C^^ is 24 eV, the [C ll] 
158 fiui ^P3/2 to ^Pi/2 line can be for med over a substantial 
part o f a H II region and according to lAannestad &c Emervl 
l|2003l ) may contribute up to 1/3 of the total intensity in the 
line. 

More recent work has indicated that approximately 25% 
of the observ ed [C ll] 158 ^m emission may come from ion- 
ized regions llAbell l 20061 ) . This is particularly true where 



low-density H ll r egions are adjacent to PDRsJ Hcilcs"l994; 
lAbel et al.i l2005l : iKaufman etal] l2006l : lAbel et al.l |20o3) . 

This effect can hamper the use of [C ll] emission as a pure 
PDR diagnostic in cases where ionized and PDR emission 
are observed in a single spectrum. Therefore, the contribu- 
tion of [C ll] emission from the ionized gas must be esti- 
mated. Such an estimate requires a separate model of the 
H II region, although in recent years computional methods 
exist that allow the H ll re gion and PDR spectrum to be cal- 
culat ed self-consistentelv (lAbel et al. 20051: Kaufman et al.l 
120061 ). A studv of S125 bv lAannestad fc Emervl (|2003D found 
that ~40% of the [C ll] and -20% of the [O l] line 
intensities come from the ionized regions. Therefore, even 
though [C ll] emission is widely observed and is usually op- 
tically thin, its dependence on the properties of the H ll 
region can diminish its use as a PDR diagnostic. The possi- 
bility that some of the [C ll] flux may come from H ll regions 
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Figure 5. Observed [C Il]i5g /[O ijea versus [O l]63/[0 l]i45 ratios, with uncertainties, compared to models of differing metallicity. The 
coloured circles represent models with different values of metallicity, ranging from 0.01 (light green circle), 0.25 (claret circle) and 0.5 
(turquoise circle) up to 5 times solar metallicity. 
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Figure 6. Observed [C lljiss /[O ijsa versus [O l]63/[0 l]i45 ratios, with uncertainties, compared to models of differing metallicity, with 
the C/O ratio adopted to be 0.53 for metallicities of 0.25 and above, but dropping to 0.16 as metallicities decrease to 0.03. The coloured 
circles represent models with different values of metallicity, ranging from 0.03 (light green circle), 0.25 (claret circle) and 0.5 (turquoise 
circle) up to twice solar metallicity. 



is not explicitly taken into account by the UCL_PDR mod- 
els. 

While ionized carbon, because of its ionization poten- 
tial, can be found in both neutral gas and ionized gas clouds, 
species such as ionized nitrogen, requiring an ionization 
potential of 14.53 eV, can arise only from H ll regions. Along 
with the [C ll] 158/xm emission line, the [N ll] 122-/im and 
205- /^m lines are the brightest emission lines contributing to 
the total far-infrared emission from our Galaxy (Wright et 
al. 1991). The ground state '^P term of the N"** ion is split 
into the three ^P2,i,o levels from which the 122 /.im and 
205 /xm lines arise. Therefore if a predicted H ll region value 
for the [C n]i58/[N n]i22 flux ratio is available, the observed 
[N ll] 122-/im flux could be used to indicate the amount of 



[C ll] 158-/im emission arising from H ll regions that are in 
the telescope beam. 

The COBE FIRAS Galactic Plane spectral measure- 
mon ts yielded an in tegrated value of [N n]i22/[N ll]205 = 1.5 
( Wright et al.lll99ll'l. The same [N ll] i22/fN n]205 flux ratio 
was measured bv loberst et al] (|2006l ) for the Great Carina 
nebula, which has Ne~30 cm~"^. They used photoionization 
models for the Carina nebula to estimate [C n]i58/[N n]i22 
= 1.6 for the H ll region. Table[3]lists the [N ll] 122-/xm emis- 
sion line fluxes measured by us from the ISO LWS spectra 
of 12 out of the 46 sources listed in Table [1] that showed a 
detection, together with the relative [C Il]i58/[N Il]i22 ra- 
tios and the estimated percentage of [C Il]i58 coming from 
the H II region. The 12 sources that have detections of both 
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Figure 7. Observed [C Il]i58 /[O ijea versus [O l]63/[0 l]i45 ratios, with uncertainties, compared to models of differing densities. In the 
left top panel the colored points represent models with densities in the range 10^^ n//^2xl0^ cm~^ decreasing from left to right. In 
the right top panel the colored points represent models with densities values in the range 10''^ n^f^SxlO* cm~'^, decreasing from left 
to right. In the left bottom panel the colored points represent models with densities in the range 10"^^ n^^S^gxlO^ cm~^ decreasing from 
left to right. In the right bottom panel the colored points represent models with densities values in the range 10'^ n^f s;9xl0^ cm~'^, 
decreasing from right to left. 



[C ll] 158 fim and [N ll] 122 have mean and median 
[C Il]i5s/[N Il]i22 flux ratios of 10.2 and 5.9 respectively. A 
H II region [C Il]i58/[N Il]i22 ratio of 1.6 implies that H ll 
regions contribute only 16% (mean case) and 27% (median 
case) of the overall [C ll] 158/im flux that is observed. 

We used the above predicted H ll region 
[C Il]i58/[N Il]i22 ratio of 1.6 along with the observed 
[N ll] 122-/im fluxes, to correct the observed [C ll] 158/im 
flux of these 12 sources for H ll region contributions. In 
Fig [9] we overplot our model results, with the same physical 
parameters used in Fig |4l on the observations. In the 
right-hand panel the observations have been corrected for 
the predicted H ll region contribution to the 158-/im fluxes 
using the measured [N ll] 122- /xm fluxes. As shown in Fig|9l 
the model results appears to fit the corrected observations 
better, compared to Fig |4j because the observations have 
been shifted to lower values of [C ll]i58/[0 ijea- The remain- 



ing discrepancy is attributed to [O l] 63-/im self-absorption, 
as discussed in Section 6, below. 



6 FITTING THE OBSERVATIONS 

The clear discrepancy between the location of the observa- 
tions and our models in Figures 3-8 can at least be partly 
attributed to self absor ption in the [ l] 6 3-Atm line. This 
is consistent with what iLiseau et al] (|2006l ) found: namely 
that ^^65% of observed [O l]63/[0 l]i45 emission line ratios 
from Milky Way PDRs are lower than can be explained by 
current models, which they attributed to optical depth ef- 
fects in the [O l] 63-/im line. They partially attributed the 
low observed ratios to absorption by cold foreground O" in 
front of the 63-/.tm-emitting region, although other explana- 
tions such as very optically thick [O l] line emission could 
not be ruled out. 
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Figure 8. Observed [C lljisg /[O l]63 versus [O l]63/[0 l]i45 ratios, with uncertainties, compared to models with differing cosmic ray 



ionization rates. The coloured circles represents models with cosmic ray ionization rates of 
(claret circle) and f=5xl0~^^s~^ (blue circle). 



:5xlO- 



(green circle), f =5x10 s 



Table 4. Line flux measurements, in units of 10 W m for extragalactic sources with ISO-LWS detections of the [N ll] 122-/xm 
emission line, together with the observed [C n]i58/[N Il]i22 flux ratio. Column 1 is the index number from Table [T] 



Index 


Source 


TDT No. of 
Observation 


[N ii] 122/im 


[C n]i,58 
[N n]i22 


[C II] 158 
from H II 


3 


NGC 253 


56901708 


1.72±0.18 


2.7±0.4 


-60% 


5 


Maffei 2 


85800682 


0.188±0.014 


5.7±0.7 


-28% 


6 


NGC 1068 


60500401 


0.356±0.026 


5.3±0.6 


-30% 


10 


NGC 1614 


85501010 


0.021±0.004 


9.7±2.5 


-17% 


11 


NGC 2146 


67900165 


0.152±0.006 


16.3±6.4 


-10% 


13 


M 82 


65800611 


2.20±0.17 


5.9±0.6 


-27% 


18 


NGC 4039/9 


25301107 


0.0468±0.0049 


8.08±1.09 


-20% 


26 


Cen A 


63400464 


0.170±0.019 


16.2±2.4 


-10% 


27 


NW Cen A 


45400151 


0.097±0.011 


28.8±5.7 


-5.5% 


28 


M 51 


35100651 


0.228±0.027 


4.2±0.7 


-38% 


29 


M 83 


64200513 


0.231±0.066 


5.2±1.7 


-31% 


30 


Circinus 


10401133 


0.213±0.019 


12.2±1.5 


-13% 



The overall conclusion of iLiseau et all l|2006l ) was that 
because this ratio depends sensitively on models, [O l] emis- 
sion has a limited diagnostic value. In fact most PDR mod- 
els calculate the emergent flux from a plane-parallel slab of 
gas illuminated from one side. A galaxy has many PDRs at 
all orientations, and optical depth effects are non- negligible. 
In the approximation that the densest PDRs correspond to 
the shells of molecular clouds and that the [C ll] line emis- 
sion and far-IR dust continuum emission are optically thin, 
while the [O l] 63-^m line is optically thick, then the [O l] 
63-fj,m line will be seen only from the front side of each cloud 
while the [C ll] line arises from both the front and the near 
sides. The velocity dispersion from cloud to cloud, however. 



allows most [O l] 63-^m photons that have escaped their 
parent cloud to escape the galaxy entirely. This scenario 
implies that we should observe only a fraction of the [O l] 
flux and all of the [C ll] flux expected from PDR models 
iMalhotra et al.ll200ll '). 

The observed [C n]i58/[0 ijes and [O l]63/[0 1)145 ratios, 
in the selected range, do not appear to be traced properly 
by our UCL_PDR model results although there is a common 
trend in the behaviour of these ratios, indicating that the ob- 
servations are off-set compared to our model results. At first 
glance the model in the lower right panel of Fig 2] where we 
varied the radiation field strength, appears to show the best 
fit to the observations. However, taking into account the ef- 
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Figure 9. Observed [C Ii]i58 /[O ijea versus [O i]63/[0 i]i45 ratios, with uncertainties. In the left panel the observations do not take 
into account the likely percentage of [C ll] coming from H ll regions, while in the right panel the estimated contribution to the [C ll] 
flux is taken into account for each source. The coloured points represent represent models with radiation field strengths in the range 
10^Go^2xlo4 and the arrows indicate the direction of the increase of radiation field strength. The arrow on the right-hand side of the 
right panel corresponds to the median shift in the [C n]i58/[N n]i22 ratio, corresponding to 27% of the [C ll] emission originating from 
H II regions (see Section 5). 



fects of [O l] self-absorption and the [C ll] contribution from 
H II regions, the models in the lower left panel of Fig 4 could 
also probably provide a fit to the observations. In fact both 
these effects would shift the models downwards and to the 
right relative to the observations. A similar trend could be 
found when we varied the density. Although at first glance 
the models in the top right panel of Fig [7] do not appear to 
show a good fit to the observations, if one takes into account 
the effects of [O l]63 self-absorption, which will shift observa- 
tions downwards and to the left relative to the models, due 
to the decrease in the [O i] 63-/im flux, these models will 
probably provide a better fit to the observations. The mod- 
els will be shifted downwards and to the right relative to the 
observations. Unfortunately the results obtained varying the 
cosmic ray ionization rate and the metallicity suggest that 
the [C ll]i58/[0 l]63 versus [O l]63/[0 l]i45 ratio is not a good 
tracer for these parameters. We identified the best parame- 
ter values of our source sample as 10''^n_r/^;9xl0^ cm~^ and 
60<Gos;8xlO^ A„= 10. These values were chosen while tak- 
ing into account the effects of absorption in the [O l] 63- /xm 
line (See Section 6.1) and the contribution to [C ll] 158-/im 
from H II regions (Section 5). 



6.1 [O l] 63- /^m self-absorption 

The intensity of the [O l] 63/xm line can be suppresse d 
through self-absorption (e.g. iGonzalez- Alfonso et al.ll2003 ). 
The "^Pi state, the upper level of the [O i] 63-//m line, is 
228 K above the g round state and its critical d ensity is ricrit 
= 4.7x10^^ cm-3 (|Tielens fc HoUenbachl [iQsBI ) . Because of 
this, at typical interstellar cloud temperatures all oxygen 
atoms occupy the lowest level of their ground state, with 
J=2. This means that [O l] 63-/xm emission, originating from 
collisional excitation in a warm medium, could be absorbed 
by the abundant reservoir of ground state neutral oxygen 



that exists within the galaxy along the line of sight. The 63- 
/im line can also be absorbed by cold foreground material 
but this is less likely in the case of a face-on galaxy, where 
the neutral oxygen column density will be less along the 
line of sight than for an edge-on galaxy. For our sample of 
sources the majority are not edge-on galaxies, therefore we 
attribute most of the suppression of the [O l] 63- /xm line to 
self-absorption in the PDRs in which the emission is excited. 

6.1.1 The [O ij 63-nm line profile 

In an attempt to investigate the effect of [O l] 63-/im 
self-absorption on the observations we have developed a 
new version of the Sph erical Muffi-Mol code (SMMOL, 
iRawlings fc Yated HooJ)). Originally SMMOL only consid- 
ered a sphere of material illuminated by the standard ISRF 
at the outer boundary of the sphere: we have now imple- 
mented the ability to insert a central illuminating source 
with an arbitrary SED. SMMOL uses an accelerated A- 
iteration (ALI) method to solve multilevel non-LTE radia- 
tive transfer problems of gas inflow and outflow. The code 
computes the total radiation field and the level populations 
self-consistently. At each radial point, SMMOL generates 
the level populations, the line source functions and the emer- 
gent spectrum from the cloud surface. This can then be 
convolved with the appropriate telescope beam. A detailed 
description of the SMM OL radiative trans f er mo del can be 
found in the appendix of lRawlings fc Yated (|200ll 'l . The cou- 
pling between the UCL_PDR code and the radiation trans- 
fer code has been performed through an interface that 
will be presented in a forthcoming paper (Bayet et al. 2010, 

^ See website: https: //www. astro .uni-koeln.de /projects 
/schilke/ sites/www. astro. uni -koeln.de. projects. schilke/ files/ 
Viti_Cologne09.pdf 
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in prep). We used the physical parameters from one of the 
UCL_PDR modelsEI and we chose an arbitrary distance of 
3.2 Mpc for the clouds modeled. In Fig[lO]we show the SM- 
MOL [O l] emission line profiles for a PDR with a diameter 
of 20 pc (left panel) and a PDR with a diameter of 5 pc (right 
panel). In Fig [TO] the effect of self-absorption is clearly seen 
in the [O iJesMm profile (solid line) while there is no trace 
of self-absorption in the [O l]i45;jm line profile (dotted line). 
The [O ijespm emission line is still significantly brighter than 
the [O l]i45,jm emission line, consistent with current obser- 
vations. The line profile from a single 20 pc GMC exhibits 
extremely strong [O l] 63-/im self-absorption. Without self- 
absorption in the [O l] 63-/im line, the 63/145 flux ratio 
would be 46.2. With self-absorption, the flux ratio is found 
to be 8.5, a reduction of 0.73 dex. The effect is smaller for 
the 5 pc model; the 63/145 flux ratio is reduced from 12.6 
to 7.9, i.e. by 0.20 dex. The log(63/145) - 0.9 flux ratios 
predicted by both SMMOL models are in reasonable agree- 
ment with the observed ratios in e.g. Fig |9l Although our 
20 pc and 5 pc diameter SMMOL models produce similar 
final 63- and 145-/im fluxes, their relative masses (1.0x10^ 
and 1.6x10'^ Mq, respectively) would seem to make an en- 
semble of the smaller clouds more plausible for matching the 
observed ratios. 

Note that although all profiles in FigllOlhave been mod- 
eled with the same physical parameters, the difference be- 
tween the two [O ijesfiin emission line profiles is due to the 
difference between the inner and outer temperatures of the 
two clouds. The 20 pc GMC is approximately 100 K cooler 
than the 5 pc GMC at the centre. This is because the 20 pc 
GMC extends to A„=19 mag, while the 5 pc GMC has a 
visual extinction of ~5 mags. These differences account for 
the much stronger [O l]63fim self-absorption in the left panel. 
The [O l]i45fim line, on the other hand, is approximately the 
same for both GMCs. 

In order to qualitatively understand the behaviour of [O 
i]63Mm self-absorption in PDRs, we investigated its sensitiv- 
ity to variations in the density and radiation field strength. 
We find that an order of magnitude increase in the radiation 
field strength causes the inner and outer temperatures to be 
twice as high as those in the models plotted in Fig. 10. This 
produces a higher flux and an increase in self-absorption 
in the [O iJesMm line- Increasing the density by one order 
of magnitude causes an increase in the optical depth of the 
cloud, so that the degree of [O ijesfim self-absorption is larger 
than in the right-hand panel of Fig llOl while the overall flux 
is reduced by an order of magnitude. We underline here that 
we are not attempting to model the actual line profiles that 
may be exhibited by these galaxies: in order to do so one 
would need to consider density gradients, multiple sources 
and the geometry of the galaxy. Our theoretical line profiles 
simply demonstrate the effect of the initial conditions on the 
degree of self-absorption in the [O ijesfiin line. 

The [O l] 63/145 flux ratios predicted by our SMMOL 
models shift them downwards by 0.2-0.7 dex relative to the 
UCL_PDR predictions, bringing them into quite good agree- 
ment with the observed ratios plotted in Figs. 9 and 10. 
However, allowance for the 0.2-0.7 dex reduction in the [O l] 



Model result listed in Fig|4](bottom left panel): f =5x10 s ^ , 
10^ cm-3, Z/Zq=1, 200Go and t=10^ yrs. 



63- /xm line strength caused by self-absorption would also 
move the UCL_PDR [C ii]/[0 ijespm model ratios in those 
figures to larger values. We ran SMMOL for the case of C^ 
and we found, as expected, that the [C ll]i58Aim emission 
line profile did not show self-absorption in any of the mod- 
els. In summary, we find that 20-80% of the intensity of the 
[O l]e3)im line can be suppressed through absorption, with 
the percentage depending strongly on the physical parame- 
ters of the PDR region. Comparing with the model trends 
shown in Fig. 9, it would seem that a combination of high 
radiation field strengths and allowance for [O l] 63-jj.m self- 
absorption could provide the best match to the observed line 
ratios. [O l] line profile observations from the Herschel Space 
Observatory should allow a greater degree of discrimination 
amongst possible models. 



7 CONCLUSIONS 

To fit the observed [C ll] 158 ^tm, [O l] 63 /.tm and [O l] 
145 ^m emission line fluxes for 28 extragalactic sources, 
measured from archival ISO-LWS spectra, we used a grid 
of 1702 PDR models from the UCL_PDR code spanning a 
large range of densities, radiation field strengths, metallici- 
ties and cosmic ray ionization rates. We took into account 
the contribution to the observed [C ll] fluxes from H ll re- 
gions, using expected [C n]i58/[N Il]i22 ratios for H ll regions 
to correct the [C lljisg fluxes of the 12 sources that had [N ll] 
122 fim detections. 

We found that the best fitting PDR models had 
lO^sCnH^gxlO" c m~^ and 60s;Go^8x l0^. Consistent with 
the conclusions of iLiseau et al.l (|2006l ). we found that the 
persistent discrepancy between the observed and modelled 
line flux ratios can at least partly be attributed to self- 
absorption in the [O l] 63-/im line. We used the SMMOL 
code to predict oxygen emission line profiles for several PDR 
models and found clearly self-absorbed [O l] 63-fim profiles, 
with 20-80% of the intensity suppressed, depending on the 
physical parameters of the PDR regions. A combination of 
high radiation field strengths and [O l] self-absorption would 
appear to provide the best explanation for the observed [C ll] 
158 /im, [O l] 63 /im and [O l] 145 line fiux ratios. 



ACKNOWLEDGEMENTS 

We thank the referee, Dr N. Abel, for constructive comments 
that helped improve the paper. MV would like to thank Dr. 
Estelle Bayet for her help throughout this research. Samuel 
Farrens is thanked for comments. 



REFERENCES 

Aannestad P. A., Emery R. J., 2003, A&A, 406, 155 
Abel N. P., 2006, MNRAS, 368, 1949 

Abel N. P., Ferland G. J., Shaw G., van Hoof P. A. M., 

2005, ApJ, 161, 65 
Abel N. P., Sarma A. P., Troland T. H., Ferland G. J., 

2007, ApJ, 662, 1024 
Asplund M., Grevesse N., Sauval A. J., AUende Prieto C, 

Blomme R., 2005, A&A, 431, 693 



Testing PDR models against ISO fine structure line data for extragalactic sources. 



13 




Velocity (km s 



Velocity (km s ) 



Figure 10. Theoretical oxygen line profiles produced from a 20 pc GMC (left panel) and from a 5pc GMC (right panel), observed at a 
distance of 3.2 Mpc. The dotted line represents the [O l]i45nm emission line profile and the solid line represents the [O i]63m™ emission 
line profile. The physical parameters of the UCL_PDR model used for modeling the oxygen line profiles were (J=5xl0"l''s"\ 10^ cm , 
Z/Z0=1, 200Go and t=10'^ yrs. 



Bayet E., Gerin M., Phillips T. G., Contursi A., 2004, 
A&A, 427, 45 

Bell T. A., Hartquist T. W., Viti S., Williams D. A., 2006, 

A&A, 459, 805 
Bell T. A., Roueff E., Viti S., Williams D. A., 2006, 

MNRAS, 371, 1865 
Bell T. A., Viti S., Williams D. A., Crawford I. A., Price 

R. J., 2005, MNRAS, 357, 961 
Bergvall N., Masegosa J., Ostlin G., Cernicharo J., 2000, 
A&A, 359, 41 

Brauher J. R., Dale D. A., Helou G., 2008, ApJS, 178, 280 
Clegg P. E., Ade P. A. R., Armand C., Baluteau J.-P., Bar- 
low M. J., Buckley M. A., Berges J.-C., 1996, A&A, 315, 
L38 

Colbert J. W., Malkan M. A., Clegg P. E., Cox P., Fischer 
J., Lord S. D., Luhman M., Satyapal S., Smith H. A., 
Spinoglio L., Stacey G., Unger S. J., 1999, ApJ, 511, 721 

Contursi A., Kaufman M. J., Helou G., HoUenbach D. J., 
Brauher J., Stacey G. J., Dale D. A., Malhotra S., Rubio 
M., Rubin R. H., Lord S. D., 2002, ApJ, 124, 751 

Crawford M. K., Genzel R., Townes C. H., Watson D. M., 
1985, ApJ, 291, 755 

Draine B. T., 1978, ApJ, 36, 595 

Fischer J., Shier L. M., Luhman M. L., Satyapal S., Smith 
H. A., Stacey G. J., Unger S. J., 1996, A&A, 315, L97 

Gonzalez- Alfonso E., Smith H. A., Fischer J., Cernicharo 
J., 2004, ApJ, 613, 247 

Habing H. J., 1968, BuU. Astron. Inst. Netherlands. , 19, 
421 

Heiles C, 1994, ApJ, 436, 720 

Higdon S. J. U., Devost D., Higdon J. L., 
Brandl B. R., Houck J. R., Hall P., 2004, 
The Astronomical Society of the Pacific, 116, 975 

HoUenbach D. J., Tielens A. G. G. M., 1997, ARA&A, 35, 
179 

Israel F. P., White G. J., Baas F., 1995, A&A, 302, 343 
Izotov Y. I., Thuan T. X., 1999, ApJ, 511, 639 



Kaufman M. J., Wolfire M. G., HoUenbach D. J., 2006, 
ApJ, 644, 283 

Liseau R., Justtanont K., Tielens A. G. G. M., 2006, 

A&A, 446, 561 
Luhman M. L., Satyapal S., Fischer J., Wolfire M. G., 

Sturm E., Dudley C. C, Lutz D., Genzel R., 2003, 

ApJ, 594, 758 

Malhotra S., Kaufman M. J., HoUenbach D., Helou G., 
Rubin R. H., Brauher J., Dale D., Lu N. Y., Lord S., 
Stacey G., Contursi A., Hunter D. A., Dinerstein H., 2001, 
ApJ, 561, 766 

Mao R. Q., Henkel C, Schulz A., Zielinsky M., Mauers- 

berger R., Storzer H., Wilson T. L., Gensheimer P., 2000, 
A&A, 358, 433 
Martm S., Martm-Pintado J., Mauersberger R., 2006, 

A&A, 450, L13 
McCaU B. J., Huneycutt A. J., Saykally R. J., Geballe 

T. R., Djuric N., Dunn G. H., 2003, The Nature, 422, 500 
Meijerink R., Spaans M., Israel F. P., 2006, ApJ, 650, L103 
Negishi T., Onaka T., Chan K.-W., Roellig T. L., 2001, 

A&A, 375, 566 
Oberst T. E., Parshley S. C, Stacey G. J., Nikola T., Lohr 

A., Harnett J. I., TothiU N. F. H., Lane A. P., Stark A. A., 

2006, ApJ, 652, L125 
Rawlings J. M. C, Yates J. A., 2001, MNRAS, 326, 1423 
RoUig M., Abel N. P., BeU T., Bensch F., Black J., Ferland 

G. J., Jonkheid B., Kamp I., Kaufman M. J., Le Bourlot 

J., Le Petit F., 2007, A&A, 467, 187 
Rubin R. H., Simpson J. P., Haas M. R., Erickson E. F., 

1991, ApJ, 374, 564 
Schilke P., Carlstrom J. E., Keene J., Phillips T. G., 1993, 

ApJ, 417, L67 
Seaquist E. R., Prayer D. T., 2000, ApJ, 540, 765 
Sternberg A., Dalgarno A., 1995, ApJS, 99, 565 
Tielens A. G. G. M., HoUenbach D., 1985, ApJ, 291, 747 
Unger S. J., Clegg P. E., Stacey G. J., Cox P., Fischer J., 

Greenhouse M., Lord S. D., Luhman M. L., Satyapal S., 



14 M. Vasta et al 



Smith H. A., Spinogiio L., Wolfire M., 2000, A&A, 355, 
885 

Wild W., Harris A. I., Eckart A., Genzel R., Graf U. U., 
Jackson J. M., Russell A. P. G., Stutzki J., 1992, 
A&A, 265, 447 

Wolfire M. G., HoUenbach D., Tielens A. G. G. M., 1989, 
ApJ, 344, 770 

Wolfire M. G., Tielens A. G. G. M., HoUenbach D., 1990, 
ApJ, 358, 116 

Wright E. L., Mather J. C., Bennett C. L., Cheng E. S., 
Shafer R. A., Fixsen D. J., Eplee Jr. R. E., Isaacman R. B., 
Read S. M., 1991, ApJ, 381, 200 



This paper has been typeset from a T^j^K/ ffl^^^X file prepared 
by the author. 



